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Abstract: An approach for solving scattering problems, based on two quantum field theory 
methods, the heat kernel method and the scattering spectral method, is constructed. This 
approach converts a method of calculating heat kernels into a method of solving scattering 
problems. This allows us to establish a method of scattering problems from a method of 
heat kernels. As an application, we construct an approach for solving scattering problems 
based on the covariant perturbation theory of heat-kernel expansions. In order to apply the 
heat-kernel method to scattering problems, we first calculate the off-diagonal heat-kernel 
expansion in the frame of the covariant perturbation theory. Moreover, as an alternative 
application of the relation between heat kernels and partial-wave phase shifts presented 
in this paper, we give an example of how to calculate a global heat kernel from a known 
scattering phase shift. 
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1 Introduction 

In this paper, based on two quantum field theory methods, heat-kernel method [1] and 
scattering spectral method [2], we present a new approach to solve scattering problems. This 
approach is a series of approaches for scatterings rather than a single approach. Concretely, 
our key result is an explicit relation between partial-wave scattering phase shifts and heat 
kernels. By this result, each method of calculating heat kernels leads to an approach of 
calculating phase shifts; or, in other words, the approach converts a method of solving 
heat kernels into a method of solving scattering problems. Many methods for scattering 
problems can be constructed by this approach, since the heat-kernel theory is well studied 
in both mathematics and physics and there are many mature methods for the calculation 
of heat kernels. 

Phase shift. All information of an elastic scattering process is embedded in a scattering 
phase shift. This can be seen by directly observing the asymptotic solution of the radial wave 
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equation. For spherically symmetric cases, the asymptotic solution of the free radial wave 
equation, £ (r 2 Jy) + Ri = k 2 Ri, is Ri (r) ' =°° (1/kr) sin (kr - Itt/2) and the 

asymptotic solution of the radial wave equation with a potential, — X Jp (r 2 jy) + + V ( r ) 

k 2 Ri, is 


Ri 


D / \ r^-oc 1 

Ri (r) = — sin 
kr 


kr-^+St (k) 


( 1 . 1 ) 


This defines the partial-wave phase shift 5i (k), which is the only effect on the radial wave 
function at asymptotic distances [3]. Therefore, all we need to do in solving a scattering 
problem is to solve the phase shift Si (k). 

Heat kernel. The information embedded in an operator D can be extracted from a heat 
kernel K (f; r, r') which is the Green function of the initial-value problem of the heat-type 
equation ( dt + D) <f> = 0, determined by [1] 


(dt + D) K (t; r, r') = 0, with K (0; r, r') = S (r — r') . (1.2) 


The global heat kernel K (t) is the trace of the local heat kernel K (t; r, r'): K (t) = 
f drK (t; r, r) = i e~ Xnlt , where A n ; is the eigenvalue of the operator D. 

The main aim of the present paper is to seek a relation between the partial-wave phase 
shift Si (k) and the heat kernel K (t; r, r 7 ). By this relation, we can explicitly express a 
partial-wave phase shift by a given heat kernel. There are many studies on the approximate 
calculation of heat kernels [1, 4-13] and each approximate method of heat kernels gives us 
an approximate method for calculating partial-wave phase shifts. 

The present work is based on our preceding work given in Ref. [14], which reveals a 
relation between two quantum field theory methods, the heat-kernel method [1] and the 
scattering spectral method [2], In Ref. [14], using the relation between spectral counting 
functions and heat kernels given by Ref. [15] and the relation between phase shifts and 
state densities given by Ref. [2], we provide a relation between the global heat kernel and 
the total scattering phase shift (the total scattering phase shift is the summation of all 
partial-wave phase shifts, S (k) = ^ (2 1 + 1) Si (k)). 

Nevertheless, the result given by Ref. [14] — the relation between total scattering 
phase shifts and heat kernels — can hardly be applied to scattering problems, since the 
total phase shift has no clear physical meaning. 

To apply the heat-kernel method to scattering problems, we in fact need a relation 
between partial-wave phase shifts (rather than total phase shifts) and heat kernels. In the 
present paper, we find such a relation. This relation allows us to express a partial-wave 
phase shift by a known heat kernel. Then all physical quantities of a scattering process, 
such as scattering amplitudes and cross sections, can be expressed by a heat kernel. 

To find the relation between partial-wave phase shifts and heat kernels, we will first 
prove a relation between heat kernels and partial-wave heat kernels. The heat kernel 
K (t; r, r') is the Green function of initial-value problem of the heat equation (1.2) with 
the operator D = —V 2 + V (r) and the partial-wave heat kernel J\) (f; r, r') is the Green 
function of initial-value problem of the heat equation (1.2) with the radial operator Di = 
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— (r 2 ^) + + V (r). By this relation, we can calculate a partial-wave heat kernel 

Ki (f; r, r') from a heat kernel K (f; r, r') directly. 

The main aim of this paper is to explicitly express the partial-wave phase shift by a 
given heat kernel. As mentioned above, by our result, each method of calculating heat 
kernels can be converted to a method of calculating scattering problems. 

In order to calculate a scattering phase shift from a heat kernel, we need off-diagonal 
heat kernels (i.e., heat kernels). For this purpose, in the following, we first calculate an off- 
diagonal heat-kernel expansion in the frame of the covariant perturbation theory. It should 
be pointed out that many methods on the calculation of diagonal heat-kernel expansions in 
literature can be directly apply to the calculation of off-diagonal heat kernels. 

A method for calculating scattering phase shifts based on the covariant perturbation 
theory in the heat-kernel theory is established as an example of our approach. 

Furthermore, we compare the scattering method established in this paper, which is 
based on the covariant perturbation theory of heat kernels, with the Born approximation. 
The comparison shows that the scattering method based on the covariant perturbation 
theory is a better approximation than the Born approximation. 

Besides applying the heat-kernel method to scattering problems, on the other hand, 
by the method suggested in the present paper, we can also apply the scattering method to 
the heat-kernel theory. In this paper, we provide a simple example for illustrating how to 
calculate a heat kernel from a known scattering result; more details on this subject will be 
given in a subsequent work. The value of developing such a method, for example, is that 
though it is relatively easy to obtain a high-energy expansion of heat kernels, it is difficult 
to obtain a low-energy heat-kernel expansion. With the help of scattering theory, we can 
calculate a low-energy heat-kernel expansion from a low-energy scattering theory. 

The starting point of this work, as mentioned above, is a relation between the heat- 
kernel method and the scattering spectral method in quantum field theory. The heat-kernel 
method is important both in physics and mathematics. In physics, the heat-kernel method 
has important applications in, e.g., Euclidean field theory, gravitation theory, and statistical 
mechanics [1, 13, 16-18]. In mathematics, the heat-kernel method is an important basis of 
the spectral geometry [1, 19]. There is much research on the calculation of heat kernels. 
Besides exact solutions, there are many systematic studies on the asymptotic expansion of 
heat kernels [20], such as the covariant perturbation theory [10-13]. With various heat- 
kernel expansion techniques, one can obtain many approximate solutions of heat kernels. 
Scattering spectral method is an important quantum theory method which can be used to 
solve a variety of problems in quantum field theory, e.g., to characterize the spectrum of 
energy eigenstates in a potential background [2] and to solve the Casimir energy [21-25]. 
The method particular focuses on the property of the quantum vacuum. 

In Sect. 2, we find a relation between partial-wave phase shifts and heat kernels. 
As a key step, we give a relation between partial-wave heat kernels and heat kernels. In 
Sect. 3, based on the relation between partial-wave phase shifts and heat kernels given in 
Sect. 2, we establish an approach for the calculation of partial-wave phase shifts, based 
on an heat-kernel expansion, the covariant perturbation theory. In Sect. 4, a comparison 
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of the approach established in the present paper and the Born approximation is given; in 
particular, we compare these two methods through an exactly solvable potential. In Sect. 5, 
we give an example for calculating a heat kernel from a given phase shift. Conclusions and 
outlook are given in Sect. 6 . Moreover, an integral formula and two integral representations 
are given in Appendices A and B. 

2 Relation between partial-wave phase shift and heat kernel: calculating 
scattering phase shift from heat kernel 

The main result of the present paper is the following theorem which reveals a relation 
between partial-wave scattering phase shifts and heat kernels. This relation allows us to 
obtain a partial-wave phase shift from a known heat kernel directly. By this relation, what 
we can obtain is not only one method for scattering problems. It is in fact a series of 
methods for scattering problems: each heat-kernel method leads to a method for solving 
scattering problems. 

Theorem 1 The relation between the partial-wave scattering phase shift, Si ( k ), and the 
heat kernel, K (■ t ; r, r 7 ) = K (f; r, 8, ip, r 7 , 9' , p'), is 

/*oo i pc~\~ioo 1 

5i(k) = 2ir 2 r 2 dr d cos yP; (cos 7 ) — / dt—K s (f; r, 9, p, r, 9’, ip') + Si ( 0 ), 

Jo J -1 Jc-io o J 

( 2 . 1 ) 

where K s (t; r, r 7 ) is the scattering part of a heat kernel, Pi (cos 7 ) is the Legendre polynomial, 
and 7 is the angle between r and r' with cos 7 = cos 9 cos 9' + sin 9 sin 9' cos (<p' — ip). 

Notice that only the radial diagonal heat kernel, K s (t; r, 6, tp, r, O', p'), appears in Eq. 
(2.1). The heat kernel K (t; r, r 7 ) is split into three parts: K (t; r, r 7 ) = I\ s (t; r, r 7 ) + 
K b (t]Y,r') + RJ (t; r, r 7 ). The free part of a heat kernel K? (t; r, r 7 ) = (t; r, r 7 ) is 

the heat kernel of the operator D = —V 2 ; the bound part of a heat kernel corresponds 
to the bound state, if exists, of the system, which, in the spectral representation, is 
K b (t; r, r 7 ) = Abound states e ( r ) V’a ( r 0; ^ Ie scattering part of a heat kernel corre¬ 

sponds to the scattering state of the system, which, in the spectral representation, is 
A s (f;r,r 7 ) = Escaping states (r) ^ (r 7 ) [14]. Note that S t (0) = n/2 if there is a 

half-bound state and Si (0) = 0 if there is no half-bound state [14]. 

The remaining task of this section is to prove this theorem. In order to prove the 
theorem, we need to first find a relation between partial-wave heat kernels and heat kernels. 

2.1 Relation between partial-wave heat kernel and heat kernel 

As mentioned above, the heat kernel K (t; r, r 7 ) of an operator D is determined by the heat 
equation (1.2) [1]. For a spherically symmetric operator D, the heat kernel K (t; r, r 7 ) = 
K (t;r,6,(p,r',9',ip') can be expressed as 

K (t; r, r 7 ) = ^ e" A " 7 Vn/m (r) V nlm ( r 0 , 

n,l,m 
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where A n i and ( r ) = Rnl ( r ) Ylm (0, T) are the eigenvalue and eigenfunction of D , de¬ 
termined by the eigenequation Di/j n i m = Xnl'^nlm- where R n i (r) is the radial wave function 
and Yf m (6,<p) is the spherical harmonics. The global heat kernel is the trace of the local 
heat kernel K (f; r, r 7 ): 


K(t) 


drK (t: r, r) = Y^ e Xnlt . 

n,Z 


(2.3) 


The local partial-wave heat kernel 


Ki (t; r, r') = Y e~ x - lt R nl (r) R nl (r 7 ) 

n 


(2.4) 


of the operator D is the heat kernel of the 1-th partial-wave radial operator [14] 

d 


}_d_ 
r 2 dr 


D ‘ = -^^Tr 


l (l + 1 ) . . 

+ o + V (r) 


(2.5) 


which determines the radial equation D[R n i = A n lRnl- The global partial-wave heat kernel 
is the trace of the local partial-wave heat kernel J\/ ( t ; r, r'), 


Ki 



r 2 drKi (t; r, r) = e 

n 


( 2 . 6 ) 


Now we prove that the relation between Ki (f; r, r 7 ) and K (t: r, r 7 ) can be expressed as 
follows. 


Lemma 2 The relation between the partial-wave heat kernel Ki ( t ; r, r') and the heat kernel 
K (i; r, r 7 ) = K (i; r, 9 , <£>, r 7 ,0 7 , <£> 7 ) is 

I\i (t]r,r') = 2 it J dcos'yPi (COS 7 ) K (t; r, 9, ip, r\ 6', ip'^j (2-7) 

and 

K (t; r, 0, ip, r', 9', ip') = Y ( 21 + !) p l (cosy) 7Q (t; r, r 7 ) . (2.8) 

Proof. In a scattering with a spherically symmetric potential, the scattering wave function 
ipnim ( r ) V 5 ) = Rnl i r ) p im (0, V 5 )- Then, by Eq. (2.2), the heat kernel can be expressed as 

i 

I< (t; r, 0 , ¥>, r 7 , 0 7 , ^ (r) (r 7 ) Y ¥>) ¥>') • (2-9) 

n,l m=—l 

Using the relation [26] 

1 2i _i_ i 

( 0 , y>) ( 0 7 , ip') = —— Pi (cos 0 cos 0 7 + sin 0 sin 0 7 cos (<p 7 — ip)) , ( 2 . 10 ) 

m=—l 

we obtain 

77 (t]r,9,ip,r',9',(p') = ^ (21 + 1) P/ (cos 7 ) y~] e~ Klt R n i (r) R ni (r 7 ) . (2.11) 

l n 
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Then, by Eq. (2.4), we prove the relation (2.8). 

Multiplying both sides of Eq. (2.8) by Py (cosy) and then integrating cosy from — 1 
to 1 give 


J d cos 'yPy (cos 7 ) K ( t ; r, 9, <p, r',9 ', ip') 


s£< 2!+1 > 


d cos 7 Py (cos 7 ) Pi (cos 7 ) 


i-i 


Ki(t;r,r'). (2.12) 


Using the orthogonality of the Legendre polynomials [26] 

f 1 2 

J d cos iPy (cos 7 ) Pi (cos 7 ) = — — - 

we obtain 




f 1 1 

J d cos 7 P// (cos 7 ) K (t; r, 6, <p, r’, 9’, cp') = —Ky (t; r, r') . 
This proves the relation (2.7). ■ 


(2.13) 


(2.14) 


2.2 Proof of Theorem 1 

Now, with Lemma 2, we can prove Theorem 1. 

Proof. In Ref. [14], we prove a relation between total phase shifts and global heat kernels, 

1 rc-\-ioo Ty~s (j.\ 

S(k) = — - { -^e k t dt + 5(0), (2.15) 

^ J c—zoo ^ 

and a relation between partial-wave phase shifts and partial-wave global heat kernels, 

1 rc-\-ioo jy-s / j.\ 

Si (k) = - / + 5, ( 0 ). (2.16) 

Here the global heat kernel and the global partial-wave heat kernel are split into the scat¬ 
tering part, the bound part, and the free part: K (t) = K s (t) + K b (t) + K? ( t ) and 
K l {t) = Kf(t) + Kf(t) + K[ (t) [14]. 

Starting from the global partial-wave heat kernel given by Eq. (2.6) and using the 
relation between partial-wave heat kernels and heat kernels given by Lemma 2, Eq. (2.7), 
we have 

/♦OO 

K t (t) = / r 2 drK 1 (t; r, r) 

Jo 

/*oo 

= 2-jt J r 2 dr J d cos 7 P; (cosy) K (t; r, 9, ip, r, 9', <p') . (2-17) 

Substituting Eq. (2.17) into Eq. (2.16) proves Theorem 1. ■ 

It should be noted here that the relations given by Ref. [14], Eqs. (2.15) and (2.16), only 
allow one to calculate the total phase shift <5 (k) from a heat kernel K (t) or to calculate the 
partial-wave phase shift 6i (k) from a partial-wave heat kernel Ki (t). Such results, however, 
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are not useful in scattering problems, because the total phase shift <5 (k) is not physically 
meaningful and the partial-wave heat kernel Ki (t) is often difficult to obtain. 

Nevertheless, the result given by Theorem 1, Eq. (2.1), allows one to calculate the 
partial-wave phase shift 5i ( k ) from a heat kernel K (t) rather than a partial-wave heat 
kernel Ki (t). The heat kernel has been fully studied and there are many well-known results 
[!]• 


3 Heat-kernel approach for phase shift: covariant perturbation theory 

In this section, based on the heat-kernel expansion given by the covariant perturbation 
theory [10-12], by the relation between partial-wave phase shifts and heat kernels given by 
Eq. (2.1), we establish an expansion for scattering phase shifts. The covariant perturbation 
theory is suitable for our purposes, since it provides a uniformly convergent expansion of 
heat kernels [10, 27]. 

The covariant perturbation theory type expansion for a partial-wave scattering phase 
shift is 5i (k) = 6 ^ (k) + (k) + • • • with 

roc 

( fe ) = - ^ J o rdrV ( r ) J h-i/2 ( kr ), ( 3 - 1 ) 

_2 poo pr 

S\ 2) (k) = -- j rdrJi +l/2 (kr) Y l+1/2 (kr) V (r) j r'dr'J? +1/2 (kr') V (r') , (3.2) 

where J v (z) and Y u (z) are the Bessel functions of the first and second kinds, respectively 
[26]. 

A detailed calculation is as follows. 


3.1 Covariant perturbation theory for heat-kernel expansion 

The heat-kernel expansion is systematically studied in the covariant perturbation theory 
[10-12], The heat-kernel expansion given by covariant perturbation theory reads [4, 6 ] 

K (t; r, r') = K < 0) (i; r, r') + it « (f; r, r') + K^ (t; r, r') + • • • 

roo 

= (r| e~ Hot + (— t) / da\da 2 5 (1 — a\ — ot 2 ) e ~ aiHot Ve~ a2Hot 

Jo 

/•OO 

+ (- 1 ) 2 / daida 2 da 3 d (1 - on - a 2 - a 3 ) e~ aiHot V e ~ a2Hot V e ~ a3Hot + ■ • • |r'>, 

Jo 

(3.3) 


where 


K {0) (t; r,r') = (r| e~ Hot |r'> = 


1 


0 -|r-r'| 2 /(4t) 


(3.4) 


(4vrt ) 3/2 

is the zero-order (free) heat kernel. Substituting the zero-order heat kernel (3.4) into Eq. 
(3.3), we obtain the first two orders of a heat kernel, 


(t\ r, r') = (r| (— t) / da\da 2 5 (1 — aq — 02 ) e aiHot Ve a2Hot \Y) 

Jo 

i 2 N 


= - / dr d 3 yV( y) 


exp 


3yrl r -y| 


4 (t-r) 


2 ) exp 


-37iy- r 


„/|2 


[ 47 T (t — r )] 3 ^ 2 


(47TT 


.3/2 


(3.5) 
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and 


(t; r, r 7 ) = (r| (— t) 2 / da\da2da^5 (1 — a\ — «2 — 03 ) e aiHot Ve a2Hot Ve a3Hot |r 7 ) 

Jo 


= J d 3 yV (y) j d 3 zV (z) J dr j dr' 

1 


exp 


-j^lr-yl 2 ) exp 


4 (t—t') 


y - z| j exp ( -jp |z - r 


„/12 


(3.6) 


[ 47 r (t — r )] 3 / 2 [ 47 T (r — t 7 )] 3,/2 ( 47 tt 7 ) 3//2 

For the spherical potentials V (r) = V (r), K W (t;r, r 7 ) and A'( 2 ) (t; r, r 7 ) given by Eqs. 
(3.5) and (3.6) become 

/*oo /* rt 

K (1) (t; r, r', 7 ) = - / y 2 ch/E (y) / / dr 

Jo J Jo 

exp (- 4 ^ (r 2 + ?/ 2 - 2 ry cos 7ry )) exp (__L ( r ' 2 + y 2 - 2 r'y cos 7r , y )) 


[4-7T (f — t)] 3 ^ 2 


(47TT 


,3/2 


(3.7) 


and 


A' (2) (t; r, r 7 , 7 ) = / y 2 A/F (y) / dQ y / z 2 cAF (2) 
jo J Jo 

rt n ' exp ^— 4 ^- 3 ^ (r 2 + y 2 — 2ry cos 7ry ) ^ 


X / dr dr' 

J Jo Jo 


(3.8) 


[4-7T (f — t)] 3//2 

ex P (^4(r-r') + " 2 - ‘ 2 V Z cos 7y z)) e xp (-^ ( z 2 + r 72 - 2zr 7 COS 7zr /) ) 


[47T (r — t 7 )] 3 / 2 


(47TT 7 


,3/2 


(3.9) 


where 7 is the angle between r and r 7 , 7ry is the angle between r and y, 7r / y is the angle 
between r 7 and y, 7yz is the angle between y and z, and 7zr / is the angle between z and r 7 . 


3.2 First-order phase shift dj * ( k ) 

In this section, we calculate the first-order phase shift in the frame of the covariant pertur¬ 
bation theory. 

The first-order phase shift <5/ (k) can be obtained by substituting the first-order heat 

kernel given by covariant perturbation theory, Eq. (3.7), into the relation between partial- 
wave phase shifts and heat kernels, Eq. (2.1), and taking r 7 = r (radial diagonal): 


^ (1) ( k ) = I r z dr-^~. I dt— I dr I y l dyV (y) 


1 [ r 2 dr— [ 

Jo 27 ™ Jc 


c+ioc e k 2 t pt 


t 


f 


exp 


r 2 W 

4p-r) 


exp 


r 2 +y 2 
4 T 


[47T (t — r)] 3//2 (47rr) 3 ^ 2 

(3.10) 


where X\ is an integral with respect to the angle, 

ry 


T\ = j d cos 7 P; (cos 7 ) J dn y exp 


r y 

2 (f — r) oos 7 r y J exp cos 7r / y 


(3.11) 
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To calculate T \. we use the expansion [26] 

OO /- 

giz cos a = £ ( 2/ + 1) i l J^-J l+1/2 (z) P t (COS a) 
Z=0 V 

to rewrite Xi as 


(3.12) 


Zi =0 


12=0 


x J d cos 7 +) (cos 7 ) J dn y P h (cos 7 ry)P Z2 (cos 7 r /y) , (3.13) 

where jL (z) = \JnJj2z) J v +\/2 ( 2 ) is the spherical Bessel function of the first kind. Without 
loss of generality, we choose r 7 = (r 7 , 0, 0) and then 7 r / y = 6 y and 7 = (9 r . Now, the integral 
with respect to Q y can be worked out directly by using the integral formula (A.l) given in 
A: 

J dQ y P h (cos 7 r y ) Pi 2 (cos 7 r / y ) = J dQ y P h (cos 7 ry ) P h (cos0 y ) = P h (cos 9 r ) ^ h,h- 

(3.14) 

The integral with respect to 7 (= 9 r ). then, can also be worked out by using the orthogo¬ 
nality of the Legendre polynomials dxPi> (x) Pi (x) =2/ (2 V + 1) [26]: 


d cos 9 r Pi (cos 9 r ) Pi x (cos 9 r 


Air 


i-i 


2 /i + l 




8tt 


(21 + 1 ) 




(3.15) 


By Eqs. (3.14) and (3.15), we achieve 


Zi = £ (2/i + 1) i ll j h 

h=o 

= 8 t ri 2l ji(- ry 


ry 


i2 (t — t) 

a 


) E ( 2! = + !) (Sr) 


j2(t — t) J \i2tJ 
Substituting Eq. (3.16) into Eq. (3.10) gives 


-I rc+ioo k 2 t ft foe 

5f’(k) = -2ir 2 — / dt—— dr y 2 dyV (y)l 2 , 


where 


To = 87 ri 


rc 

i 21 / 
Jo 


r 2 dr 


exp 


27ri 


4(i-r) 


exp 


4t 


/o [47T (t — r )] 3//2 (47rr ) 3//2 \i2(t t) J \i2r) 

To calculate the integral X 2 , we use the integral representation, Eq. (B.l), given in B to 
represent the factor ji ( i2 [ t l T ) ) + {+?+) as 


ry 


Ji 


(ry\ 


(3.16) 


(3.17) 


(3.18) 


Ji 


ry 


Jl 


( ry\ = 1 


sin 


i2(t — t) J \i2r J 


dcosd- 


+ (Jf) 2 - 2 ^L^cos 0 

i 2 (t— t) \i 2 t / i 2 (t—T) i 2 t 


1-1 


' 2 + (3f) 2 -^Wt+,+r^0 


Pi (cos 0 ). 


ry 


i 2 (t—r) 


(3.19) 



























Substituting the integral representation (3.19) into Eq. (3.18) and working out the integral 
give 


r*l 


Z -2 = 47n' 2/ / d cos 9Pi (cos ( 

J -1 Jo 


fC 

< 0) 

Jo 


r 2 dr 


exp 


r 2 +y 2 
4 (t~r) 


exp 


r 2 +y 2 
4 T 


[47T (t — t)] 3 / 2 (47rr) 


x3/2 


sin 




ry 

i2 (t—r) 




r y 

i2{t—r) 


•21 


r*l 


-/ d cos 6 Pi (cos 0) exp-(cos 0 + 1)1. 

(47Tt) 3/2 J-l V / 


(3.20) 

(3.21) 


Substituting Eq. (3.21) into Eq. (3.17) and performing the integral with respect to r, we 
have 

«P> <*) - -£** [v*iW (V) 4- £2° £ ic0B „, (cos9) exp (_££ 

(3.22) 

Using the expansion exp (— y 2 cos 0/ (2f)) = ££ =0 (21 + 1) i l ji (^i £ ) Pi (cos 0) (see Eq. (3.12)) 
and the orthogonality of the Legendre polynomials, we can work out the integral: 


nl , m / y 2 cos 0 

d cos 9Pi (cos 0) exp I — 


2 1 


y£ (21' + l) i l ji/ f i— ^ f d cos 6 Pi (cos 0) Pj/ (cos 0) 
i'=o \ 2t JJ -1 

•2Z 2>\fTrt ,. 


= 

J ' 2t 


= % 


b+i/2 


2t 


(3.23) 


where I v (z) is the modified Bessel function of the first kind and the relation ( 2 ) = 
\f^TJ2z)i l I u -\-\/2 (~) is used. Substituting Eq. (3.23) into Eq. (3.22), we have 


£OC 

*l(k) = -- ydyV (y) 
z Jo 


1 /■c+ioo /„,2 

' . * V 6 "’' /(2,)/ < + l/2 

irv'i ^ 


27ri 


(3.24) 


Finally, by performing the inverse Laplace transformation in Eq. (3.24), 

1 /-c+ioo fc 2 t / 2\ 

2Si £ ioo /<2,), «/2 ( s ) - +V2 (*r). ( 3 . 25 ) 

the first-order phase shift given by covariant perturbation theory, Eq. (3.1), is obtained. 


/o') 

3.3 Second-order phase shift S { (k ) 

In this section, we calculate the second-order phase shift in the frame of the covariant 
perturbation theory. 

The second-order phase shift 5) (k) can be obtained by substituting the second-order 

heat kernel given by covariant perturbation theory, Eq. (3.9), into the relation between 
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partial-wave phase shifts and heat kernels, Eq. (2.1), and taking r' = r: 

poo i pc-\-ioo &k 2 t rt pt 

<5 , 2 (k) = 27t 2 / r 2 dr -— T / dt - / dr / dr' 

Jo 27u Jc-ioc t J 0 J 0 


POO POO 

x / y 2 dyU (y) / z 2 cPU (z) 
jo jo 


exp exp (-Jg±£y) exp (-p±F) 


[47T (t — r)] 3 ^ 2 [47T (r — r')] 3 ^ 2 (47rr') 3 ^ 2 


X 3 , 

(3.26) 


where 
X, = 


•3 = J dCOS'}Pi (cos 7 ) J dQy j cALexp 

Using Eq. (3.12), we rewrite the integral X 3 as 


exp 


yz cos 7 y2 


2 (r — r 7 ) 


TV exp 


zr cos 7 zr ' 
2 U 
(3.27) 


X 3 = ^( 2/ 1 + l)^j h 


7=0 


ry 


*2 (t — r) 


7=0 


( 2/2 + 0 *‘ 2 ^ ( i2 (t- rU ) ^ + 0 ih Jh ( 


7=0 


i2t' 


x J d cos 7 ^ (cos 7 ) J dtty J dkl z P h (cos7 ry ) P / 2 (cos7 yz ) P 7 (cos7 zr /) . (3.28) 

Without loss of generality, we choose r 7 = (r 7 , 0, 0) and then we have 7 zr / = 0 Z . The integral 
with respect to Q z can then be worked out by use of the integral formula, Eq. (A.l), given 
in Appendix A: 

J dtt z Pi 2 (cos 7 y z) Pi 3 (cos 7 zr /) = j dPl z P i2 (cos 7 yz ) Pi 3 (cos 9 Z ) = Pi 2 (cos 9 y ) 

(3.29) 

The integral with respect to £l y also can be integrated directly by Eq. (A.l), 

/ 47T 47T 47T 

dn y P h (cos 7 r y ) Pi2 (cos 6 y) + 1 ^ 2 ,i3 = Ph (cos 9 r ) + ^ 7,7 9/ + ^ 7,7- (3.30) 

Then, performing the integral with respect to 7 (7 = 0 r when r 7 = (r 7 ,0,0)) in Eq. (3.28), 
we have 

47 T 47T /* 3 327T 2 

2/i + ^ 7,7 2l2 + -^ 7,7 J ^dcos0 r Pi (cos Q r ) P h (cosfl r ) = + 3 fyA, fA,7 ■ ( 3 - 31 ) 

By Eqs. (3.29), (3.30), and (3.31), we have 

2,37; f r 2/ 


X 3 = 32TT 2 P l ji — 


J2 {t — t) ) 2 1 \i 2 (r — t 7 ) y/ Jl \i2r' 
Substituting Eq. (3.32) into Eq. (3.26) gives 

1 /■c+ioo fc 2 t rt rr roo 

<*>j 2 W = 647T 4 i 31 — / dt— dT dr' / y 2 dyV (y) 
J c —zoo t Jo Jo Jo 

.. r ^ - 7 v~ 


yz 


zr 


(3.32) 


/ 

JO 


x / r 2 dr 


r 2 +y 2 


00 2 , exp J ex P (- 17 ' 


[47T (t — r)] 3//2 (47TT') 3 / 2 \i2(f t)^ Vz2r 7 




Jz T 


zr 


(3.33) 

(3.34) 
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To perform the integral with respect to r, by using the integral representation (B.l) given 
in B, we rewrite 


Jl 


ry 


i2 {t — t) ) 7 1 \i2r') 


. f zr \ If 1 


Sill 


dcos9- 




ry 

i2(t—r) 


'-1 


ry 


i2 


i2{t— t) 


+ (^) 2 -2^^cos0 


Pi (cos 9 ). 


(3.35) 


Then, the integral with respect to r can be worked out, 


r 2 . ex P (-w£)) ex p(-^) . ( ry \ . ,zrs 

Jo r r [4tt (t - r )] 3/2 (4vr r ') 3/2 31 \i2(t - t)) 31 VttT') 


1 

2 

sin 


roo exp 
2 


r- 2 +y 2 

4(t-r) 


exp 


_zj±rj 

4r' 


/ d cos 9Pi (cos 0 ) / r 2 dr- , , 

U Jo [4vr (t - r)] 3/2 (47 tt') 3/2 




i2(t—r) 


’ 2 +(!&) 2 -2s^^c OS « 


ry 
i2(t—r) 


12 o 

I _ 2_n/_ rn sd 

+ l *2r' > Z i2(t—r) %2 t' 0 


If 1 

= — / d cos 9P/(cos 8) -- 7 - 

8W-! a ^ [ 47r(t _ r + T /)]3/2 

Substituting Eq. (3.37) into Eq. (3.34), we have 


exp 


y 2 + z 2 + 2 yz cos 9 
4 (t — T + t') 


(3.36) 

(3.37) 


. . 1 /*C+200 t /* 

5 2) (k) = 8n 3 i 31 - [ dt - [ 

1 27xiJ c _ ioo t J 0 


diP t I dr I dr ' / y 2 dyV (y) 


0 Jo 


l 


r 00 exp 

x / z 2 dzV (z) 


4(r—r') J 

7 Jl 


0 

exp 


[47 t (r — r')] 3//2 (t t 7 ) 


yz 


.V 2 +z 2 


4(i-r+r'); f 1 ( yz COS 0 

- 77777 / a cos tfn cos d ) exp- 7 -— 

[ 47r ( t _ r + r /)] 3 /2 y _! U ' P V 2 ( t-r + r ') 


(3.38) 


(3.39) 


Using the expansion (3.12) and the orthogonality of the Legendre polynomials, we have 

yz cos 9 


J d cos 9 Pi (cos 9) exp 

= y, i. 21 ' + !) iV Ji 


2 (t — T + T 1 ) 


v =0 


= 2 i l ji I -- 


yz 


i 2 (t - t + t') J J_i 


d cos 9Pi (cos 9) Pi / (cos 9) 


yz 


i2 (t — t + t') J 

Substituting Eq. (3.41) into Eq. (3.39) and setting T = r — t' , we have 

1 rc+ioo k 2 t roo roo 

4 ■ (k) = 16tt 3 — / dt— / y 2 dyV (y) / z 2 dzV (z) 

zm Jc-ioo t Jo Jo 


(3.40) 

(3.41) 


rt rr exp 

x dr dT 
Jo Jo 


y 2 +z 2 


■ - 1 ,1 y 2 +z 2 \ 

4 T ) CX P \ 4(t-T)) ■ ( yz \ ■ f V z 


(4vr Tf 2 [47r(t-T)f 2Jl ^TJ Jl { i2 (t — T) 


(3.42) 
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Exchanging the order of integrals Jjj dr J Q r dT —> dT ftp dr and resetting T = t' give 


rc+ioo gk 2 t fOC 


5^ (fe) = 167 T 3 -—: / dt - / y 2 dyV (y) / z 2 dzV (z) 


2iri 


dr / dr- 


y 2 +z 2 


, ft exp exp 


;; 2 +z 2 


(47TT') 3 / 2 [47T (t — T 7 )]" 

Integrating with respect to r, we have 


yz \ ■ 


_ r /')l 3 / 2 ^ \ i2 (f — r 7 ) 


yz 


(3.43) 


where 


c(2) 


-1 /*00 /*00 1 /*C+200 

<fc) =s 7 0 v‘ d « v wj ^vw—j 


z dzV ( 2 ) x—: / die Z 4 , 

C—200 


(3.44) 


I 4 = - I dr - 


exp 


y 2 +z 2 
4(i—r') 


exp 


y 2 +^ 2 

4r' 


(* - T ') 


1/2 


r'3/2 




yz 


i 2 (t — t') 


Ji 


yz 

i2t' 


= -8 


, ,, -fc 2 t f k2 ji (%) (%) ff ( kz ), v > z 
Jo X k2 j 2 (ky )31 (kz) n i (kz ), y < z ’ 


(3.45) 


where n; (z) is the spherical Bessel function of the second kind. Thus, the inverse Laplace 
transformation of Z 4 can be worked out: 


1 /»C+200 

2 i 7 ri J q — 200 


dte k ~ t l4 = —4 


k 2 ji ( k y) ni (ky) jf (kz ), y > z 
k 2 j 2 (ky) jl (kz) ni (kz ), y < z 


Substituting Eq. (3.46) into Eq. (3.44), we have 


6 i 2) (k) = ~k 2 I y 2 dyji (ky) n t (ky) V (y) I z 2 dzj'f (kz) V (z) 


r 00 

u2 / n 2 


'0 

r»oo 


(3.46) 


/*oo /*oo 

fc 2 / y 2 dyjf (fey) 14 (y) / £ 2 dzj; (fez) ?i; (kz) V (z) . (3.47) 

do Jv 


By exchanging the order of integrals, J 0 °° dy f°° dz —>• J 0 °° dz J 0 ~ dy, we rewrite Eq. (3.47) 


as 


(fe) = -A : 2 / y 2 dyji (ky) n t (ky) V (y) I z 2 dzj'f (kz) V (z) 

Jo Jo 

r»oo 


roo 

u 2 / ,. 2 . 


rv 


roo rz 

k 2 / z 2 dzj; (fez) u; (fez) 14 (z) / y 2 dyj 2 (ky) V (y ). (3.48) 

do do 


Obviously, the two parts in Eq. (3.48) are equal. Using j; (z) = y 7 r/ (2z)J ; +1 / 2 (z) and 
71; (z) = \/ 7 r/ ( 2 z)y; + i /2 (z) gives Eq. (3.2). 


4 Comparison with Born approximation 

The approach for scattering problems established in the present paper is to convert a method 
of calculating heat kernels into a method of solving scattering problems. As an application, 
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in Sect. 3, we suggest a method for the scattering phase shift, based on the covariant 
perturbation theory of heat kernels. 

In scattering theory, there are many approximation methods, such as the Born ap¬ 
proximation, the WKB method, the eikonal approximation, and the variational method 
[ 281 . 

In this section, we compare our method with the Born approximation. 


4.1 Comparison of first-order contribution 

For clarity, we list the result given by the above section in the following. 

The first-order phase shift given by covariant perturbation theory given in Sect. 3 reads 

roo 

6 l ( k )cpt = ~2 J 0 rdr ' J l+l/2 ( kr ) V ( r ) ■ ( 4J ) 

For comparison, the first-order phase shift given by the Born approximation reads [28] 

■ POO ' POO 

d l (k) B om = arctan -- / rdrV (r) jf +1/2 (kr) ~ -- / rdr jf +1/2 (kr) V (r) + • • • . 

(4.2) 

Obviously, the leading contributions of these two methods are the 
proximation, the first-order contribution is in fact arctan — (7r/2) J 0 °° 
but the higher contribution can be safelv ignored in the first-order cc 


same (in the Born ap 


die me &dine ^111 tiie uuin cx± 

L v-/ 2 ) / 0 °° rdrV ( r ) J l+ 1/2 ( kr ) 
irst-order contribution). 


^lUAimaiiun, niiDi-uiuu iS in fact arctan 

but the higher contribution can be safely ignored in the 

4.2 Comparison of second-order contribution 

The second-order phase shift given by covariant perturbation theory given in Sect. 3 [Eqs. 
(3.2) and (3.48)] reads 

2 roo rr 

<r Wept = - y J rdrJ l+1/2 (kr) Y l+1/2 (kr) V (r) j r'dr'Jf +l/2 {kr') V (F) 

2 roc rr 

-j ^ rdrJ i+ i/a (kr) Y l+1/2 (kr) V (r) J ^ r'dr'jf +1/2 {kr) V (/) 

2 poo poo 


ir 2 r°° 


JO 

POO 

rdrJ i+i /2 ( kr ) V ( r ) / r ' d r'J l+ 1/2 {kr') Y l+1/2 ( kr') V {r') . (4.3) 

^ JO Jr 

The second-order phase shift given by the Born approximation [28] reads 

r 2 poo pr 

6 i 2) (k )Born = arctan rdrJ l+1/2 (kr) Y l+1/2 (kr) V (r) J ^ r'dr'Jf +1/2 {kr') V {r') 

2 roc roo 

~ J o rdrJ i+ 1/2 ( kr ) v (r) J r'dr' J l+1/2 {kr') Y l+1/2 {kr') V {r') 

2 roc rr 

- J 0 rdr J i +\/2 {kr) Y l+1/2 {kr) V (r) j ^ r'dr'J? +1/2 {kr') V (r') 

~2 roo roo 

- — y rdrJ 2 +1/2 (kr) V (r) J r'dr'J l+1/2 {kr') Y l+l/2 {kr') V {/) + ■■■ . 

(4.4) 

It can be directly seen that the leading contribution of the second-order Born approxi¬ 
mation and the leading contribution of the second-order covariant perturbation theory are 
the same. 
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4.3 Comparison through an exactly solvable potentialczU (r) = a/r 2 

In this section, we compare the two methods, the covariant perturbation theory method, 
and the Born approximation, through an exactly solvable potential: 

V(r) = ± (4.5) 

Using these two approximation methods to calculate an exactly solvable potential can help 
us to compare them intuitively. 

The phase shift for the potential (4.5) can be solved exactly, 


Si 




(4.6) 


In order to compare the methods term by term, we expand the exact result (4.6) as 5i = 
<5^ + sj 2 ^ + • • •, where 


(1) = ira 

1 2 ( 2 /+ 1 )’ 

(2) = 7ra 2 

1 2 ( 2 /+ 1 ) 3 ' 


(4.7) 

(4.8) 


First order: The first-order contribution given by covariant perturbation theory and 
the Born approximation can be directly obtained by substituting the potential (4.5) into 
Eqs. (3.1) and (4.2), respectively: 


4° m„„, = - 


ira 


c Pt 2(2/+ 1)’ 


<*, (1) (k) 


Born 


= arctan 


7ra 

ira 1 

ira 

2(2/+ 1) 

2(2/ + 1) 3 

2(2/+ 1) 


(4.9) 

(4.10) 


Comparing with the direct expansion of the exact solution, Eqs. (4.7) and (4.8), we 
can see that both results are good approximations, and the result given by covariant per¬ 
turbation theory is better than the result given by the Born approximation. 

Second order: The second-order contribution given by covariant perturbation theory 
and the Born approximation can be directly obtained by substituting the potential (4.5) 
into Eqs. (3.2) and (4.4), respectively: 


+ Mot = 


7i a 


c p* 2(2/ + l) 3 ’ 


Si 2) (k) 


Born 


= arctan 


ira 2 

7ra 2 1 

7ra 2 

2(2/+ l) 3 

2(2/ + l) 3 3 

2(2/ + l) 3 _ 


(4.11) 

(4.12) 


Comparing with the second-order contribution, Eq. (4.8), we can see that, like that 
in the case of first-order contributions, the result given by covariant perturbation theory is 
better. 
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5 Calculating global heat kernel from phase shift 


The key result of this paper is a relation between partial-wave phase shifts and heat kernels. 
Besides solving a scattering problem from a known heat kernel, obviously, we can also 
calculate a heat kernel from a known phase shift. Here, we only give a simple example 
with the potential a/r 2 . A systematic discussion of how to calculate heat kernels and 
other spectral functions, such as one-loop effective actions, vacuum energies, and spectral 
counting functions, from a solved scattering problem will be given elsewhere. 

For the potential 

rv 

(5.1) 


V{r) = ^ 


the exact partial-wave phase shift is given by Eq. (4.6), 


Si = ~ 


7r 


l + 


+ OL — I l + 


(5.2) 


By the relation between a global heat kernel and a scattering phase shift given by Ref. 


[14], 


2 r°° 

Kf (t) = -t / kdkS t (k) e 
n Jo 


—krt 


Si( o) 


7r 


we can calculate the scattering part of the global heat kernel immediately, 


*W = -2 


' a+[i+ D ~{ i+ \ 


(5.3) 


(5.4) 


In this case, the bound part of heat kernel Kf (t) = 0 and the free part of heat kernel 
Kj (t) = R/y/Airt — 5 (l + 5 ): where R is the radius of the system. The global partial-wave 
heat kernel then reads 


K l (t) = Knt) + K*(t) + K{ (t) 

R 1 I Z 1 

\Z4tt t 2 V \ 2 


(5.5) 


As a comparison, we calculate the partial-wave heat kernel for V (r) = a/r 2 by another 
approach. 

The partial-wave heat kernel of a free particle, I\(. which is the heat kernel of the 
radial operator D^ ree = — (r 2 ^) + , can be calculated directly: 


K ‘ (t;r ’ r ' ) = 2i^ eXP 


1 f r 2 + r' 2 


4 1 


1 1 + 1/2 


rr 
2 1 


(5.6) 


By setting + fz, where s = \Ja + {l + 1/2 ) 2 — 1/2, we can obtain the 

partial-wave heat kernel of operator D[ = — \ {r 2 -^J) + ^ 3 ^ + 


5 


Ki (t; r, r') = 


2 t\frr' 


exp 


r 2 _|_ r /2 


4t ) y/ a +H+ 1 / 2) 2 l 2t 




(5.7) 
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Taking trace of A) (f;r, r') gives the global partial-wave heat kernel A) (t): 


Ki (t) 


,-R 


r 2 drKi (t; r, r) 


r2(i+v) 


(4f) 1+r? T (2 + 77 ) 


2F2 



r] + 1 ; ?7 + 2 , 2 ?? + 1 ; 



(5.8) 


where 77 = \J{l + 1/2 ) 2 + a, p F q (ai, 02 • • ■ a p ; &i, 62 ■ ■ ■ b q ; 2 ) is the generalized hypergeomet¬ 
ric function [26]. Expanding A") (t) at A —> 00 gives 


(t) 


R 1 ^ 2 t 7 2 - i + ie~ R 

—== - —77 + vt—- 2 -—- 1 - 

2 1 V 4^/ttR 


(5.9) 


When A —>• 00 , one recovers the heat kernel given by Eq. (5.5). 


6 Conclusions and outlook 

In this paper, based on two quantum held theory methods, the heat kernel method [1] and 
the scattering spectral method [ 2 ], we suggest an approach for calculating the scattering 
phase shift. The method suggested in the present paper is indeed a series of different 
methods of calculating scattering phase shifts constructed from various heat kernel methods. 

The key step is to find a relation between partial-wave phase shifts and heat kernels. 
This relation allows us to express a partial-wave phase shift by a heat kernel. Then, each 
method of the calculation of heat kernels can be converted to a method of the calculation 
of phase shifts. 

As an application, we provide a method for the calculation of phase shifts based on the 
covariant perturbation theory of heat kernels. 

Furthermore, as emphasized above, by this approach, we can construct various methods 
for scattering problems with the help of various heat kernel methods. In subsequent works, 
we shall construct various scattering methods by using various heat-kernel expansions. 

In this paper, as a byproduct, we also provide an off-diagonal heat-kernel expansion 
based on the technique developed in the covariant perturbation theory for diagonal heat 
kernels, since the heat kernel method for scatterings established in the present paper is 
based on the off-diagonal heat kernel rather than the diagonal heat kernel. It should be 
emphasized that many methods for calculating diagonal heat kernels can be directly applied 
to the calculation of off-diagonal heat kernels. That is to say, the method for calculating the 
diagonal heat kernel often can also be converted to a method for calculating off-diagonal 
heat kernels and scattering phase shifts, as we have done in the present paper. Therefore, 
we can construct scattering methods from many methods of diagonal heat kernels, e.g., 
[1, 9, 29], 

The heat kernel theory is well studied in both mathematics and physics. Here, as 
examples, we list some methods on the calculation of heat kernel. In Refs. [30-32], the au¬ 
thors calculate the heat-kernel coefficient with different boundary conditions. In Ref. [33], 
using the background field method, the author calculates the fourth and fifth heat-kernel 
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coefficients. In Refs. [34-36], the authors calculate the third coefficient by the covariant 
technique. In Refs. [9, 37], by a string-inspired worldline path-integral method, the au¬ 
thors calculate the first seven heat-kernel coefficients. In Ref. [29], a direct, nonrecursive 
method for the calculation of heat kernels is presented. In Ref. [38], the first five heat-kernel 
coefficients for a general Laplace-type operator on a compact Riemannian space without 
boundary by the index-free notation are given. In Refs. [10-13, 27, 39-43], a covariant 
perturbation theory which yields a uniformly convergent expansion of heat kernels is es¬ 
tablished. In Refs. [44-46], a covariant pseudo-differential-operator method for calculating 
heat-kernel expansions in an arbitrary space dimension is given. 

An important application of the method given by this paper is to solve various spectral 
functions by a scattering method. The problem of spectral functions is an important issue 
in quantum field theory [15, 47, 48]. A subsequent work on this subject is a systematic 
discussion of calculating heat kernels, effective actions, vacuum energies, etc., from a known 
phase shift. We will show that, based on scattering methods, we can obtain some new 
heat-kernel expansions. It is known that though there are many discussions on the high 
energy heat-kernel expansion, the low-energy expansion of heat kernels is relatively difficult 
to obtain. While there are some successful low-energy scattering theories, by using the 
relation given in this paper we can directly obtain some low-energy results for heat kernels. 

Starting from the result given by the present paper, we can study many problems. 
The method presented in this paper can be applied to low-dimensional scatterings. One- 
and two-dimensional scatterings and their applications have been thoroughly studied, such 
as the transport property of low-dimensional materials [49-51]. We will also consider a 
systematic application of our method to relativistic scattering. The relativistic scattering 
is an important problem, e.g., the collision of solitons in relativistic scalar field theories 
[52] and the Dirac scattering in the problem of the electron properties of graphene [53, 54], 
We can also apply the method to low-temperature physics. There are many scattering 
problems in low-temperature physics, such as the scattering in the problem of the transition 
temperature of BEC [55, 56] and the transport property of spin-polarized fermions at low 
temperature [57, 58]. 

The application of the method to inverse scattering problems is an important subject 
of our subsequent work. The inverse scattering problem has extreme significance in physics 
[59, 60]. In practice, for example, the inverse scattering method can be applied to the 
problem of BEC [61] and the Aharonov-Bohm effect [62], 

In Ref. [47], we provide a method for solving the spectral function, such as one-loop 
effective actions, vacuum energies, and spectral counting functions in quantum field theory. 
The key idea is to construct the equations obeyed by these quantities. We show that, 
for example, the equation of the one-loop effective action is a partial integro-differential 
equation. By the relation between partial-wave phase shifts and heat kernel, we can also 
construct an equation obeyed by phase shifts. 

Moreover, in conventional scattering theory, an approximate large-distance asymptotics 
is used to seek an explicit result. In Ref. [63], we show that such an approximate treatment 
is not necessary: without the large-distance asymptotics, one can still rigorously obtain an 
explicit result. The result presented in this paper can be directly applied to the scattering 
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theory without large-distance asymptotics. 


A f dQ'Pi (cos 7 ) Pi 1 (cos 9') 

In this appendix, we provide an integral formula: 

f A/tt 

/ dtt'Pi (cos 7 ) Pu (cos O') = 21 + i Pl ( cos0 ) hi'■ (A- 1 ) 

where 7 is the angle between r = (r, 9 , <f>) and r 7 = (r 7 , 6' , <j>') and dVt' = sin O'dQ'dcf)'. 

Proof. Using the integral formula [64] 

I dn'Y l0 ( 7 ) y vo (e 1 ) = (0) hu (A. 2 ) 

and the relation Y/o (0, (j) = y 7 (2 1 + 1) / ( 47 r)P; (cos0), we have 

J dn'Yio ( 7 ) Yj/q (#') = ( cos Pv ( cos 0/ ) 

= Pi (cos 9) 5u>. (A.3) 


This proves Eq. (A.l). ■ 

B Integral representations of ji (u) ji (v) and j/ (u) rii (v) 

In this appendix, we provide two integral representations for the product of two spherical 
Bessel functions j 1 (u) ji (v) and ji (u) ni (v). 

3l (u) ji (v) = \ [ dcos6 SlI h t hp l (cos 9 ), (B.l) 

2 J -1 w 

where w = \fu 2 + v 2 — 2 uv cos 6 and l is an integer. 

Proof. Using the expansion [65] 

OO 

= V] (2Z + 1) ji (u) ji (v) Pi (cos 0) , (B.2) 

1=0 

where u = |u| and v = |v| with 9 the angle between u and v. Multiplying both sides of 
(B.2) by Py (cos0) and integrating from 0 to it give 

d cos 9 (21 + 1) ji (u ) ji (v ) Pi (cos 9) P v (cos 9) = 2 j v (u) j v ( v ). 

(B.3) 

Here, the orthogonality, J ’_ 1 d cos 9Pi (cos 9) Py (cos 9) = 2/ (21 + 1) Sw, is used. This proves 
Eq. (B.l). ■ 


[ d cos 9 5111W Py (cos 9) = 
l_ 1 w 


00 

£ 

1=0 ' 


3l i u ) n l ( v ) 


1 

2 



dcos9 


cos u; 


w 


Pi (cos 0) , 


U < V, 


(B.4) 
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where w = \Jv? + v 2 — 2uv cos 8. 
Proof. Using the expansion [65] 


cos w = _ (21 + 1) ji (u) ni (u) Pi (cos 6), u<v, (B.5) 

w 1=0 

where u = |u| and v = |v| with 8 as the angle between u and v. Multiplying both sides of 
(B.5) by Py (cos 8) and integrating from 0 to 7 r give 



d cos 8 


COS U) 


10 


Pi' (cos 8) 


oo 

£ 

(21 

+ 

1 ) ji (u) ni (v 

1=0 




£ 
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This proves Eq. (B.4). ■ 
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